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ABSTRACT 

We investigate the large-scale influence of outflows from AGNs in enriching the IGM with metals in 
a cosmological context. We combine cosmological simulations of large-scale structure formation with 
a detailed model of metal enrichment, in which outflows expand anisotropically along the direction of 
least resistance, distributing metals into the IGM. The metals carried by the outflows are generated by 
two separate stellar populations: stars located near the central AGN, and stars located in the greater 
galaxy. Using this algorithm, we performed a series of 5 simulations of the propagation of AGN- 
driven outflows in a cosmological volume of size (128ft. _1 Mpc) 3 in a ACDM universe, and analyze the 
resulting metal enrichment of the IGM. We found that the metallicity induced in the IGM is greatly 
dominated by AGNs having bolometric luminosity L > 10 9 L©, sources with I0 8 < L/Lq < I0 9 having 
a negligible contribution. Our simulations produced an average IGM metallicity of [O/H] = —5 at 
z = 5.5, which then rises gradually, and remains relatively flat at a value [O/H] = —2.8 between z = 2 
and z = 0. The ejection of metals from AGN host galaxies by AGN-driven outflows is found to enrich 
the IGM to > 10 — 20% of the observed values, the number dependent on redshift. The enriched 
IGM volume fractions are small at z > 3, then rise rapidly to the following values at z = 0: 6 — 10% 
of the volume enriched to [O/H] > -2.5, 14 - 24% volume to [O/H] > -3, and 34 - 45% volume 
to [O/H] > —4. At z > 2. there is a gradient of the induced enrichment, the metallicity decreasing 
with increasing IGM density, enriching the underdense IGM to higher metallicities, a trend more 
prominent with increasing anisotropy of the outflows. This can explain observations of metal-enriched 
low-density IGM at z ~ 3 - 4. 

Subject headings: cosmology — galaxies: active — galaxies: jets — intergalactic medium — methods: 
N-body simulations 
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1. INTRODUCTION 

Several lines of evidence, both observational and the- 
oretical, demonstrate that the accretion of matter onto 
supermassive black holes (SMBHs) at the centers of ac- 
tive galaxies, and the resulting feedback from them, have 
participated intimately in the formation and evolution of 
galaxies, and strongly influenced their present-day large- 
scale structures ( e.g.. iSafpeterl 119641; lLvnden-Belllll96£ 



Reesl 119841: iRi chstone et al.1 |1998t iFerrarese fc Merrit 
200a Kauffmann fc Haehneltjl2000T: iSazonov et al, 
Hop kins et al.ll2006|:lLapi et al.ll2006l : lMalbon et al 



200E; 



2007 



Menci et all 120081 : Tshankarl |2009D .~ A large fraction of 
AGNs are o bserved to host outflows, in a wide variety o f 
forms (e.g.. iCrenshaw et all [200llChartas et al.ll2009h : 
theoretical studies also indicate triggering of massive 
galactic winds from black holes (|King fc Pounds! 120031 : 
Monaco fe Fon tanot 2005) . Searching for AGN feedback, 



Nesvadba et al.l (|2006l) observed that the outflow from a 
powerful radio galaxy at z — 2.16 (likely a forming mas- 
sive galaxy) could remove a significant amount of gas 
(< 50%) from a L > L* galaxy within a few tens to 100 
Myr. They concluded that AGN winds might have a cos- 
mological significance comparab le to, or perhap s larger 
than, starburst-driven winds. iNesvadba et al.1 (2008) 
found spectroscopic evidence for bipolar outflows in three 
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powerful radio galaxies at z ~ 2 — 3, with kinetic ener- 
gies equivalent to 0.2% of the rest-mass of the SMBH. 
These outflows possibly indicate a significant phase in 
the evolution of the host galaxy. 

AGN outflows can enrich the intergalactic medium 
(IGM) with metals produced inside the host galaxy. It 
is of great interest to study the impact of this metal 
enrichment on the evolution of the IGM and the sub- 
sequent formation of galaxies. While the energy de- 
posited by outflows can potentially inhibit the formation 
of low-mass galaxies, metals enrichment can potentially 
counter-balance this effect by reducing the cooling rate 
of the intergalactic gas. At present, the full mechanism 
and operation of AGN feedback and metal enrichment on 
different scales is poorly understood, b oth from theoret- 
ical a nd observational points of view (jMathur fe Fields! 
l200l . 

The metallicity in the vicinity of AGNs have been ob- 
served to be super-solar, and non-evolving over a wide 
redshift range, indicated by various studies as follows. 
Observational spectroscopy of high-redshift quasars in- 
dicate that broad-line regio ns possess gas metallici- 
ties several times so l ar (e.g.. | Hamann fc Ferlandl IT9921 : 
Dietrich et all 120031 : iNagao et al.l l2006at iWang et al.l 
20091) . reaching as hi gh as ~ 15Z ([Baldwin et al.l l2003) 



Nagao et a l. (2006a) found no metallicity evolution over 



2.0 < z < 4.5 for a given quasar luminosity, and esti- 
mated the typical metallicity of br oad-line region (BLR) 
gas clouds to be Z ~ 5Z Q . IDietrich et all (2009) 
found that the gas metallicity of the BLR is ~ 3Zq. 
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Uuarez et al.l (|2009ft inferred that the metallicity of the 
BLR gas is very high (several times solar) even at z ~ 
6, and observed a lack of metallicity evolution among 
quasars within 4 < z < 6.4 and those at lower-z, as 
found in other studies. In addition to the BLRs, narrow 
absorption line systems in QSOs and Seyferts have been 
observed to have metal abundances above the solar value 
(jD'Odorico et al.ll2004t iGroves et al.ll2006[) . with the gas 
metallicity in the NLRs non-evolvin g over 1 < z < 4 
(|Nagao et al.l l2006bt iMatsuoka et~all 120091) . Cosmolog - 
ical hydrodynamic simulations (|Di Matteo et al.l 12004 ) 
also indicate that quasar hosts have supersolar (Z/Zq ~ 
2 — 3) metallicities already at z ~ 5 — 6, and the rate of 
evolution of the mean quasar metallicity as a function of 
redshift is generally fiat out to z ~ 4 — 5. 

The average metallicity of the IGM is found to rise 
above Zigm ^ 1O _3 Z by z ~ 2 — 3, as indicated by sev- 
eral studies, outlined below. Observations have detected 
carbon (C IV) lines in the low neutral hydrogen column 
density Lyman-a forest clouds toward z ~ 3 quasars. 
Clouds with TVhi > 10 15 cm -2 are fo und to have a car- 
bon abundance of ~ 10~ 3 — 10~ 2 solar (iCowie et al.lll995b 
iSongaila fc Cowid 119961: iSongailal 11997ft. Studyi ng the 
evolution of intergalactic metals. ISoneailal (|2001l ) found 
that the C IV column density distribution function is in- 
variant over the redshift range z = 1.5 — 5.5, and an IGM 
metallicity of 2 x l O~ 4 i? is already in place at z = 5. 
iRvan- Weber et all (|2009l ). surveying intergalactic met- 
als between 5.2 < z < 6.2, deduced that the CIV mass 
density at a mean (z) = 5.76 is smaller by a factor of 
~ 3.5 compared to the value at z < 4.7. Their study 
implies an IGM metallicity Zigm > 10~ 4 Z & . Oxygen 
(O VI) observations show heavy-element abundances in 
the range ~ 10~ 3 — 10~ 15 solar at z ~ 2 - 2.5 (e.g., 
Bergeron et"aT1l2002t ISimcoe et alj [2002t iCarswell et all 
2002t iSimcoe et al.l 120041 ). Comparing the results of 
numer ical simulations to ob servations, iHellsten et al.l 
(| 19971 ) and lRauch et al.l (|1997l ) showed that the observed 
column density ratios of different ionic species at z ~ 3 
can be well reproduced if a ratio [C/H] ~ —2.5 is as- 
sumed in the models, with a scatter of roughly an order 
of magnitude. Such high observed metallicities are ex- 
plained by models involving star formation early on at 
higher redshifts (z > 6), in cluding massive star formation 
(jGnedin &: Ostriken 11997ft, and SN-driven preg alactic 
outflows (|Madau et al.l 120011: iScannapieco et all 120021 ). 
Some of these earliest metals could have or iginated in 
Popu l ation III stars at very hi g h redshifts dCarr et all 
1981 IQstriker fc Gnedinl 119961 IHaiman fc Loebl 119971 : 
Wise fc Abell 120081 ). At low redshifts, O VI observations 
indicate a mean co smic metallicity of log ( Zicm /Zp)) > 
-2.4 at z = 0.9 (|Burles fc Tvtlerl 119961 ). In the lo- 
cal Universe we have [C/H] = -1.2 and [O/H] > -2 
(iTripp et al 112002ft. metallici ties less than 2.5-10% of so- 
lar values (IShull et al.ll2"0"03ft. and a mean o xygen metal- 
licity of ~ QMZ? (jDanforth fc Shullll2005l ). 

The IGM enrichment is highly inhomogeneo us. In cos- 
mological simulations, iCen fc Ostrikerl (| 19991 ) obtained 
global average metallicity values increasing from 1% of 
the solar value at z = 3 to 20% at present, and noted 
the very strong dependency on local density, with high- 
density regions having much higher metallicity than low- 
density regions. iSchave et al.l ([2003) found that the car- 



bon abundance in the IGM is spatially highly inhomo- 
geneous and is well described by a lognormal distri- 
bution for fixed overdensity and redshi ft. Regions of 
low-m etallicity have also been observed: iPrescott et al.1 
(2009) discovered a w 45 kpc Ly-a nebula at z ~ 1.67, 
exhibiting strong He II emission and very weak C IV 
and CIII] emission, which implies low metallicity gas 

(z IGM < lo- 2 - io- 3 z ). 

Several studies, going to as early as iVoitl (|1996| ). have 
showed that winds (starburst, SN-driven, or quasar) 
from early galaxi es drive shocks t hat h eat, ionize, and 
enrich the IGM. iTegmark etall (|1993ft modeled SNe- 
driven winds and predicted that most of the IGM is 
enriched to at least 10% of the current metal content 
by z = 5. Adding a prescription for chemical evolution 
and metal ejec t ion b y winds in cosmological simulations, 
lAguirre et al.l (|2001l ) calculated the enrichment of the 
IGM at z > 3 by galaxies of baryonic mass > 1O 8 ' 5 M . 
They found that winds of velocity > 200 — 300 km 
s _1 can enrich the IGM to the mean observed level, al- 
though many low-density regions wo uld remain metal- 
free. Carrying out SPH simulations, IScannapieco et ahl 
(2006) showed that the observed CIV correlation func- 
tions cannot be reproduced by models in which the IGM 
metallicity is constant or a local function of overden- 
sity. However, the observed properties are consistent 
with a model in which metals are confined within bubbles 
with a typical radius R s about sources of mass > M s , 
with best-fitting values of R s ps 2 comoving Mpc and 
M s « 1O 12 M at z = 3. 

On a different note, I Gnedinl ()1998l ) showed, using cos- 
mological simulations, that most of the metals in the 
IGM are transported from protogalaxies by the merger 
mechanism, and only a sm all fraction of metals is deliv- 
ered by SNe-driven winds. iFerrara et al.l (|2000ft showed 
that SNe explosions and subsequent blowouts fail by 
more than one order of magnitude to pollute the whole 
IGM to the metallicity levels observed, requiring some 
additional physical mechanism(s) to be more efficient 
pollutants. Employing the adiabatic feedback model, 
where stella r feedback natura lly drives winds in SPH 
simulations, iShen et aD fj2009f ) found that IGM metals 
primarily reside in the warm-hot intergalactic medium 
(WHIM) (metallicity of 0.01 - 0.1 Z Q with a slight de- 
crease at lower- z) throughout cosmic history. They also 
found that galactic winds most efficiently enrich the IGM 
fo r intermediate mass gala xies (10 10 — 1O 11 ' 5 M ). 

IScannapieco et all ()2002l ) performed Monte Carlo cos- 
mological simulations to track the time-evolution of SN- 
driven metal-enriched outflows from early galaxies, find- 
ing that up to 30% volume of the IGM is enriched to 
above 1O _3 Z at z = 3, making the enrichment biased 
to the areas near the starbursting galaxies themselves. 
The majority of enrichment occurred relatively early 
(5 < z < 12), with mass-averaged cosmological metal- 
licity values 10~ 3 — lO" 1,5 ^. U sing a hydrodynamic 
simulations, iThacker et al.1 (|2002f ) showed that outflows 
from starbursting dwarf galaxies (with total halo masses 
< 1O 1O M ) enrich ~ 20% of the simulation volume with 
a mean metallicity of 0.3% solar at z = 4 . Perform- 
ing SPH simulations, Khala tvan et al.l (|2008ft found that, 
without AGN feedback, metals are confined to vicinities 
of galaxies, underestimating the observed metallicity of 



3 



the IGM at overdensities < 10, while the AGN feedback 
model and the ef ficiency they assum e d ove renrich the 
underdense IGM. lOppenheimer et"al . (2009) performed 
high-resolution cosmological simulations to study the ob- 
servability and physical properties of five ions in absorp- 
tion between z = 8 and z = 5. The volume filling factor 
of metals they calculated increases during this epoch, 
reaching ~ 1% for Ziqm > 10~ 3 ^ by z = 5. Us ing hy- 
drodynamic simulations, iTornatore et al.l (|2010D found 
that feedback from winds driven by supernovae and BHs 
leave distinct signatures in the chemical and thermal his- 
tory of the IGM (especially at z < 3), with BH feedback 
providing a stronger and more pristine enrichment of the 
warm-hot IGM. 

While most of the previous work on dispersing met- 
als into the IGM has concentrated on winds driven by 
SNe or starbursts, some authors have focus e d on feed- 
back by AGN-driven outflows. IMoll et al.l (120071) ex- 
amined metal enrichment of the intracluster medium 
(ICM) caused by A GN outflows i n gala xy clusters, be- 
tween = 1 — 0. iFabian et al.l (|2010f) studied of the 
effect of AGN feedback on metal enrichment and ther- 
mal properties of the ICM in hydrodynamical simula- 
tions, considering feedback from gas accretion onto cen- 
tral SMBHs and that in the 'radio mode'. Several 
authors (e.g.. IGopal-Krishna fe W iita 2001; Bar ai et all 
l200i IGopal-Krishna et al.l I2004t IBarai fc Wiital 120071) 
have qualitatively argued that the huge lobes of radio 
galaxies could contribute substantially to spreading met- 
als into the IGM, by sweeping out the metal-rich ISM of 
young galaxies which they encounter while expanding. 
In a recent study, direct observational evidence for out- 
flow of m etal-enriched gas driven by a radio AGN were 
found by Kir kpatrick et all (|2009f ) . Using Chandra ob- 
servations of the Hydra A galaxy cluster, they showed 
that the metallicity of the ICM is enhanced by up to 
0.2 dex along the radio jets and lobes (enhancements ex- 
tending over a distance of 20 — 120 kpc from the central 
galaxy) compared to the metallicity of the undisturbed 
gas. They estimated that 10 — 30% of the iron mass 
within the central galaxy has been transported out. 

Oth er cosmological stud ies on quasar/radio galaxy out- 
flows (iFurlanetto fc Loeb1l200lHScannapieco fe Ohf l2004; 
iLevine fe Gnedirjl2005l: lBarail l2008) have considered out- 
flows expanding with a spherical geometry. However, in 
realistic cosmological scenarios, where the density dis- 
tribution shows significant structures in the form of fil- 
aments, pancakes, etc., outflows are expected to ex- 
pand anisotropically o n large scales. In a previous paper 
(|Germain et al.l [20091 hereafter Paper I), we designed a 
semi-analytical model for AGN outflows expanding with 
an anisotropic geometry, following the path of least resis- 
tance, which we implemented into a numerical algorithm 
for cosmological simulations. We used this algorithm to 
study the distribution and enriched volume fraction of 
metals produced by a cosmological population of AGNs, 
in a ACDM universe, over the age of the Universe. These 
results are presented in Paper I. In the second part of this 
project, we have modified the algorithm used in Paper I 
by adding a prescription for the amount and composi- 
tion of the metals transported by AGN outflows. Using 
this modified algorithm, we now complete the study of 
Paper I by calculating the metallicity and chemical com- 
position of the IGM. 



This paper is organized as follows. In we describe 
the numerical methodology for incorporating a metal en- 
richment technique into our anisotropic outflow model. 

describes the cosmological iV-body simulations we 
performed. The results are presented and discussed in 
21 We summarize and present our conclusions in fJSJ 

2. THE NUMERICAL METHOD 

2.1. The Basic Algorithm 

Our basic algorithm is described in details in Pa- 
per I. We use a Particle-Mesh (PM) code to simulate 
the growth of large-scale structure in the universe, in a 
cubic volume with periodic boundary conditions expand- 
ing with Hubble flow. At each time step, we filter the 
density distribution to identify the density peaks (local 
maxima), where AGNs will be located. We assume that 
each AGN has an active lifetime ^agn = 10 8 yr, and that 
AGN luminosities range from 10 8 L ( 7 ) to 10 14 Ly with th e 
luminosity distribution function of lHopkins et al.l {2007). 
We then determine the number of AGNs and their birth 
redshifts in order to reproduce the redshift-dependent 
AGN distribution function N(L,z). As the simulation 
proceeds, we locate, at each time step, the appropri- 
ate number of newly-born AGNs at randomly-selected 
density peaks. Using scaling relations between the mass 
Mbuige of the bulge and the mass Mbh of the central 
black hole, and between Mbh and the AGN luminosity, 
we estimate the baryonic mass M& of the host galaxy, 
and its total mass M ga i = M&f2o/^&,o> where ttb,a and 
fio are the baryon and total density parameters, respec- 
tively. We assume that a fraction /outflow = 0-6 of thes e 
AGNs produce outflows (|Gangulv fe Brothertonl 120081) . 
Our outflow mod el is based on the o utflow evolution 
equations given in lTegmark et all (I1993D. combined with 
the anisotropic outflow model of IPieri et al.l ()2007D . In 
this model outflows form two bipolar cones of opening 
angle a, expanding along the direction of least resistance 
around the density peak. We treat the opening angle 
as a free parameter. For details, we refer the reader to 
Paper I. 

We note that using an analytical outflow model nec- 
essarily requires tha t we make some simpl ifying assump- 
tions. The model of iTegmark et all (|1993f) assumes that 
the external medium has an uniform density. Following 
it, our model of outflow evolution approximates that the 
density of the environment through which the AGN out- 
flows propagate is equal to the mean baryon density of 
the universe, at the corresponding redshift. In a gen- 
eral scenario, the density of the external medium would 
vary both with distance and direction. However, the 
directional dependence of the density is taken into ac- 
count in our model, by using anisotropic outflows that 
travel along the direction of least resistance. In that 
se nse, our model c onstit utes an improvement over that 
of ITegmark et al.1 (|1993[ ). The distance dependence of 
the density is taken into account in the outfl ow expan- 
sion model used by ILevine fc GnedirJ (120051), w hich is 
also an improvement over ITegmark et al.1 (| 19931) . Our 
model does not take into account of this distance depen- 
dence. Rather, the focus of our study is on the direc- 
tional dependence of the density field, and the resulting 
anisotropic outflows. Expanding along the path of least 
resistance, our model outflows travel in regions that are 
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first overdense, then underdense, hence the effect of the 
density variations partly cancel out. 

This basic algorithm enabled us to identify the regions 
in space that were metal-enriched by outflows, and to 
calculate the enriched volume fraction (Paper I), but did 
not consider the actual amount of metals transported by 
the outflows. We have therefore modified this algorithm 
to take into account the mass and chemical composition 
of the metals contained in each outflow, in order to cal- 
culate the metal content and metallicity of the enriched 
IGM. The modifications to the basic algorithm are de- 
scribed in §§2.2-2.4 below. 

2.2. Mass of Metals in Galaxy 

Metals are generated within a galaxy by stars undergo- 
ing SNe explosions. The AGN outflows act to spread the 
produced metals to the broader IGM. We distinguish two 
regions of stellar populations as the sources of metal gen- 
eration: (1) stars near the AGN (very center of galaxy), 
and (2) stars located away from central AGN within the 
rest of the galaxy. 

The observational diagnosis of the metallicities is dif- 
ferent in the two regions, which makes it difficult to dis- 
entangle them at high redshifts. Observational studies of 
AGN/quasar metal abundances measure metallicities of 
the broad line and narrow line gas clouds near the AGN, 
which have a high metal content (region 1). Studies of 
galactic metallicity measure the abundances in the larger 
galaxy ISM (region 2), which are g enerally lower tha n in 
region 1. Previous simulations bv I Moll et al.1 (j2007j ) as- 
sumed that only the central metals (region 1) are driven 
out to the IGM by AGN outflows. They neglected the 
fact that AGN outflows can entrain some enriched ISM of 
the larger galaxy and spread it to the IGM. In this work, 
we take into account both sources of metals in regions 1 
and 2. 

Quasar abundances studies indicate that the mass of 
gas enriched in metals near the central engine is at least 
equal to the mass of t he central black hole, M gas > Mbh 
(jHamann et al.ll2007t ). In this work, we adopt the lower 
limit (M gas = Mbh) as the mass of metal-enriched gas 
near the AGN (region 1) which is carried by the outflow. 
In the absence of well-constrained values and for sim- 
plicity, we assume a central metallicity of Z\ = 5Zq for 
all the AGN in our simulations. Such an assumption 
is supported by spectroscopic observations of quasars 
(§Q}, which have detected several times solar metallic- 
ity of broad-line regions, w ith no metallicity evolution 
over a broad redshift range. iMoll et ail (|2007l ) also used 
a value of 5Z Q as their outflow metallicity. Hence, the 
mass of metals in the central region which is carried away 
by an outflow and dispersed into the IGM is given by 
Mz,x = -ZiMbh = 5ZqMbh- We take the mass frac- 
tion of heavy elements in the Sun as Z Q = 0.02 (jCoxl 
2000), to be consistent with the value used in our anal- 
ysis (detailed in Appendix A) for computing the metal 
mass fraction. Such values ar e based on older estimate s 
of the Solar metal abundance ([Anders &: Grevesselll989f ). 
and have also been used by other studies with which we 
compare our results. 

For region 2, a metallicity value is ascribed to each 
AGN host galaxy according to its mass. We adopt the es- 
ti mation of redshi f t-depe ndent mass-metallicity relation 
by iMaiolino et al.l ((20081) , where the following formula- 



tion is used, 

12+log(n /n H ) = -0.0864(logM*-logM ) 2 +^o, (1) 

where the metallicity is expressed as the logarithm of 
oxygen over hydrogen abundance (number density ratio) , 
log(n /n H ), and M* is the galaxy stellar mass0 Mq and 
Kq are fitting parameters whose values are determined at 
each redshift to obta i n best fit to the observations. Table 
5 of IMaiolino et ail ((2008) gives the parameter values 
over the redshift range z = 0.07—3.5. For a given redshift 
z, we obtain the parameters by interpolating between two 
consecutive lines in that table. For redshifts z < 0.07 
and z > 3.5, we simply use the values at z = 0.07 and 
z = 3.5, respectively. The stellar mass and ISM mass 
of the galaxy are M* = /*M 6 and Mism = (1 - /*)M b , 
respectively, where /* is the star formati o n effic iency. We 
use the value /* = 0.1, as in IPieri et al.1 ((20071) . 

Using the value of n /ri H obtained from equation ((TJ), 
we then calculate the mass ratio of metals using 

Z ° = 7 / i ' ( 2 ) 

(n /n H ) V pi . H + Mz,o 

where 

V-Z . Mpr,H = + MHc . (3 

Z>2 ° H 

and the sum in equation ((3)) is over all metals (all ele- 
ments except hydrogen and helium). The derivation of 
these relations is presented in Appendix A. 

For each AGN, the metallicity Zq of its host galaxy is 
computed using equations ((TJl-Q. Zq is then multiplied 
by the gas mass in the galaxy ISM (Mism) to get the 
mass of metals Mg j6 „ in the ISM. Only these metals 
can be carried by an outflow, since the metals locked 
up inside stars cannot escape. We assume that a fixed 
fraction / GSC of the metals in the ISM will be carried by 
the outflow (if the re is an outflow), and we use the value 
of / esc = 0.5, as in lPieri etldl ((20071) . The mass of metals 
carried away by an outflow from the non-central regions 
of its host galaxy is then given by Mz2 = foscMc.cn — 
f esc ZQ(l — f±)Mb. The total mass of metals carried by 
the outflow is therefore 

M z , out = 5Z M BH + f csc Z G {l - U)M h . (4) 

2.3. Distribution of Metals from Galaxy to IGM 

We simulate the enrichment of the IGM by distribut- 
ing the metals transported by the outflows to the parti- 
cles (of the PM code) intercepted by them. The exact 
mechanisms by which galactic metals are distributed by 
outflows to the IGM are not well-known, so we have to 
make some simplifying approximations. We assume that 
the metals carried by the outflow escape the galaxy at a 
constant rate during the active-AGN life, starting from 
the time i D ir when the outflow is born until the time t g 
when the central engine turns off. The total mass of met- 
als carried from the central AGN and the greater galaxy, 
Mz jOUt = Mz,i +Mz,2, is considered to be deposited uni- 
formly into the IGM gas overlapped by the outflow. We 
implement this into the algorithm as follows. 

4 Maiolino ct al. (2003) use the notation O/H for n /n H . 



5 



1) Each outflow will deposit its total metal content 
-Mz.out into the diffuse IGM, over the time the AGN is 
active (i.e., from birth to switch-off epoch). The time 
during which the outflow expands in the active-AGN 
phase is t AGN = t Q ff — t^ ir = 10 8 yr. The metals car- 
ried by the outflow are considered to be deposited into 
the IGM at an uniform rate in time during the active life- 
time. Hence, during the active-AGN phase, an amount 
of metals -Mz jOU t Ai/i AGN is added to the outflow during 
a timestep At. If at that time the outflow intercepts N p 
particles, the mass of metals added to each intercepted 
particle is then Mz iOU tAi/(i A Q N iVp). This insures that 
the total mass of metals deposited in the IGM during the 
active-AGN life is equal to Mz jOU t- 

2) Simply following such a prescription (point 1) leads 
to the following issue. When a metal-enriched outflow 
permeates ambient matter, the parts of IGM being per- 
meated first gain more metals than the parts permeated 
later during the expansion of the outflow. In our simu- 
lations, this adds more metals to the particles near the 
AGN (which are overlapped by the expanding outflow 
starting from a time close to tbir), than the particles lo- 
cated far away (which are intercepted later at a time 
close to i ff), as the outflow expan ds. This might indeed 
be the case in reality [the model of (|Tegmark et al.|[l993l) 
assumes a linear velocity field inside the outflow, which 
implies no mixing]. However, turbulent motion inside 
the outflow could redistribute the metals. 

In this paper, we assume that the distribution of metals 
in the outflow is uniform at all times. To implement this 
we redistribute the metals within the outflow volume, by 
performing an averaging of the metal mass contained by 
the relevant particles. At each timestep, we identify all 
the N p particles intercepted by the outflow, add up the 
metals (as explained in point 1 above), and redistribute 
these metals evenly among the N p particles. Such a pro- 
cedure smooths out the metal distribution in the IGM 
volume overlapped by the outflow, such that at every 
timestep, all particles inside the outflow have the same 
metal content, which steadily increases during the active 
AGN life as more metals are added. 

3) After the active-AGN phase is completed (at time 
t = t ff), all the metals carried by the outflow have been 
deposited into the IGM. At that point, the outflow enters 
a post- AGN dormant phase (see §2.5 of Paper I). The in- 
ternal pressure of the outflow exceeds the external pres- 
sure of the IGM, driving the expansion further. During 
that phase, more IGM gas is being swept by the outflow 
(in practice, this means more particles are overlapped 
by the outflow) , but the metal content of the outflow no 
longer increases. To implement this we follow the same 
averaging technique (detailed in point 2 above) for this 
redistribution of metals. The particles newly overlapped 
(after t g) gain a fraction of metals from the particles 
that were overlapped earlier (between t D1 r and £ ff)> such 
that the total metal mass in the outflow is conserved. 
Eventually, the outflow reaches pressure equilibrium with 
the external IGM. At that point, the outflow passively 
follows the Hubble flow, and no new IGM gas is being 
swept. 

We need to be careful in situations when outflows over- 
lap. If a particle is hit by several outflows, it should have 
a higher metal content than particles hit by only one 
outflow. However, the averaging prescription described 



above would tend to erase this, by giving to all parti- 
cles inside a same outflow the same metal content. The 
proper way to deal with this situation would be to keep 
track of the amount of metals each particle has received 
from each outflow. For instance, suppose that a particle 
has been enriched by two outflows A and B, which have 
deposited an amount of metals Mz,A an d Mz,b onto that 
particle, respectively. When averaging the metal content 
of particles inside outflow A, only the mass Mz,a would 
be involved in the averaging. Similarly, when averaging 
the metal content of particles inside outflow B, only the 
mass Mz,b would be involved in the averaging. This 
would conserve the total mass of metals, while insuring 
that particles hit by several outflows have a higher metal 
content. 

In practice, this is extremely difficult to implement in 
the algorithm. We have considered several approaches, 
and all of them require either the memory or the CPU 
time to scale like the product (number of particles) x 
(number of outflows), that is, 256 3 x 1.5 x 10 6 . To deal 
with this situation, we make the approximation that 
Mz,a = Mz,b- More generally, if a particle P over- 
laps with N outflows, we assume that each outflow con- 
tributed a fraction 1 /N of the metals contained in that 
particle. Then, when applying the averaging technique 
to one of these outflows, only a fraction 1 /N of the metals 
contained in particle P is included in the averaging. The 
result is that, within each outflow, the particles hit by 
that outflow only will have the same metal content, while 
particles hit by N outflows will contain roughly N times 
more metals. This might result in some local fluctua- 
tions when particles are overlapped by outflows having 
very different metal content, but with so many outflows 
in our computational volume, the overall effect of this 
approximation should be small. Also, a more luminous 
outflow (higher Mbh and Mb) would carry more metals 
(high Afz,out), but at the same time would grow larger, 
intercepting more particles (larger N p ) within its volume 
than a less luminous outflow. So the amount of met- 
als deposited to a particle Mz, out At /(t^^Np) remains 
comparable for outflows with different luminosities. No- 
tice that the total mass of metals is strictly conserved in 
this procedure. 

The methodology described in § §2.21 and 12.31 is imple- 
mented in the PM code, hence the propagation of out- 
flows and the enrichment of the IGM is simulated along 
with the growth of large-scale structure. In the next sec- 
tion f E]2.4[) . we describe calculations that use the dumps 
produced by the simulations (which contain the posi- 
tions, velocities, masses, and metal content of the PM 
particles). These calculations are performed as post- 
processing after the simulations are completed. 

2.4. Computation of the Metallicity of the IGM 

The outflowing matter consists of a mixture of primor- 
dial gas (H, He) and metals carried from the host galaxy. 
The mass of H and He transferred by an outflow to the 
overlapped IGM volume is assumed to be negligible com- 
pared to the mass of H and He already present in the 
IGM. Only the mass of metals deposited into the IGM is 
considered to be significant and counted. Any IGM vol- 
ume is considered to have no metals before it has been 
encompassed by an outflow. Hence, in the initial con- 
ditions, the metal content of each PM particle is set to 
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TABLE 1 

Parameters of the Simulations 



Run a (°) Lmin/i© £max/£© 



A 


180 


10 8 


10 14 


B 


120 


10 8 


10 14 


C 


60 


10 s 


10 14 


CI 


60 


10" 


10 14 


C2 


60 


10 8 


10 9 



zero. 

In order to compute the metallicity of the IGM at a 
given time (or redshift), we take the corresponding dump, 
which contains the particle positions, masses, and metal 
content. We then divide the computational volume into 
Ng — 256 3 cubic cells and we compute the matter density 
p C cii and metal density /? me tai in the center of each cell, 
using a smoothing method, as described in Paper I. We 
then multiply these numbers by the cell volume to get the 
total mass M ce n and total metal mass Mz in each cell. 
Assuming a universal ratio of baryons to dark matter, the 
mass of baryonic gas (comprising of H, He, and metals) 
in the cell is (-\o/^o)A/ceii- 

Once the AGN outflows are born and start propagating 
into the IGM, the IGM is enriched by metals carried by 
the outflows. If Mz is the mass of metals added to a 
grid cell (from outflows hitting it), the metal abundance 
ratio, by mass, of the cell is 

M z M z 
coU M H + M He + M z ~ M H + M Hc 
Mz = M Z 

~ (M H + M H e) initial (Q b , /Qo)M ceU K1 

Here it has been assumed that the mass of metals im- 
parted to a cell is negligible compared to the mass of 
H and He already there (Mz <C Mh + A/hc), and that 
the total mass of H and He in the cell remains the same 
before and after outflow impact, the added mass being 
negligible. 

Using equation ([2]), we can calculate the number den- 
sity ratio of O over H, 

W/cell Vz.O V 1 - Z cell) ' 

This gives us the oxygen abundance ratio in each grid 
cell of the computational volume. 

We finally express the metal content in terms of the 
metallicity, which is the logarithmic ratio of the metal 
abundance of a cell over the corresponding solar value: 
[0/H] ce ii = log 10 (n o /ri H ) C ei] - log 10 (n o /n H ) Q . We com- 
pute the metallicity using oxygen, as it is the most abun- 
dant element (after hydrogen and helium) in the uni- 
verse. The metallicity for any other element can be sim- 
ilarly computed, assuming some abundance ratio. 

3. THE COSMOLOGICAL SIMULATIONS 

As in Paper I, we simulate the growth of large-scale 
structure and the propagation of outflows in a cubic cos- 
mological volume of comoving size Lbox = 128 h^ 1 Mpc 
= 182.6 Mpc with periodic boun dary conditions, using 
a Particle-Mesh (PM) algorithm (|Hocknev fc Eastwood! 



11988ft . with 256 3 equal mass particles and a 512 3 grid. 
This corresponds to a particle mass m par t = 1.38 x 
10 10 M Q , and a grid spacing A = 0.357 Mpc. 

We consider a ACDM model with a present baryon 
density parameter flb,o = 0.0462, total matter (baryons 
+ dark matter) density parameter Slo = 0.279, cosmo- 
logical constant ,7a, o = 0.721, Hubble constant H Q = 
70.1kms~ 1 Mpc" 1 (h = 0.701), primordial tilt n s = 
0.960, and CMB temperature Tqmb = 2.725, consis- 
tent with the results of WMAP5 combined with the data 
from baryonic acoustic oscillations and supernova stud- 
ics__l (Hinshaw ct al. 2008) . We generate initial conditions 
at redshift z — 24, and evolve the cosmological volume 
up to a final redshift z = 0. The simulations presented 
in this paper use the same initial conditions as simula- 
tions A, B, and C of Paper I. 

We performed a series of 5 simulations by varying some 
parameters of our AGN evolution model. Table [T] sum- 
marizes the characteristics of each run. The first column 
gives the letter identifying the run. In columns 2, 3, 
and 4, we list the opening angle of the outflows, and the 
minimum and maximum bolometric luminosities used for 
generating the AGN population from the QLF (§2.2, Pa- 
per I), respectively. These five simulations were all pre- 
sented (among others) in Paper I. We redid them with 
our modified algorithm which incorporates the explicit 
metal enrichment prescription described in <j2] 

In addition, we performed an extra simulation run in- 
vestigating a different prescription for the initial location 
of the AGN population, whose results are given in 34.71 

4. RESULTS AND DISCUSSION 

The results are presented and discussed in the follow- 
ing subsections. We compute various statistics of the 
metallicity of oxygen (following the prescriptions in §2.4(1 
to quantify the resulting IGM metal enrichment in our 
simulations. 

4.1. The Final Radius of the Outflows 

The final radius R of an outflow depends on the mass 
Mg a i and birth redshift Zbir of the galaxy producing that 
outflow. Examining these dependencies will help us to 
understand the results presented in the following sec- 
tions. 

We consider Run C. Figure [T] shows the final radius of 
the outflows vs. birth redshift, for galaxies in the mass 
range M ga i = [1.0,1.1] x 1O 7 M . The thickness of the 
curve is caused by the use of a finite mass interval. The 
radius is nearly independent of redshift for outflows with 
Zbi r > 2. This may seem surprizing at first. Outflows 
that form earlier are facing a much stronger external 
IGM pressure, and thus reach pressure equilibrium much 
sooner. However, once they have reached that equilib- 
rium, they do not stop. Instead, they expand with Hub- 
ble flow. For outflows that reach equilibrium, the final 
physical radius is equal to the comoving radius at equi- 
librium. When outflows reach equilibrium early, their 
physical radius at that time is small, but their comoving 
radius is large. 

There is a slight break in the slope at redshift z = 1.1, 
followed by a rapid decline down to R = at z = 0. 

5 http: / /lambda. gsfc.nasa.gov/product/map/dr3/parameters_summary 
.cfm 
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Fig. 3. — Quantities in the whole (182.6 Mpc) 3 simulation vol- 
ume as a function of redshift. The top panel shows the total bary- 
onic mass of all host galaxies (Y2 Mb, dot-dashed upper curve), and 
the total mass of metals transported by outflows from all galaxies 
(V) Mz out, solid curve). The total metal mass is a summation 
of metals carried from the central regions (Y2Mz.i, dotted curve), 
and metals carried from the non-central regions (the greater ISM) 
(Y] Mz.2, dashed curve). The middle panel gives the ratio of metal 
mass over galaxy mass. The upper two panels are the same for runs 
A, B and C. The bottom panel shows the fractional volume of the 
simulation box enriched by the outflows in run C (a = 60°). 

ever, this dependence is very weak. We find that R scales 

1/5 

approximately like M ^ . This is a result of competing 
effects. Following equation (9) of Paper I, the thermal 
pressure generated by the AGN luminosity is what drives 
the expansion of the outflow. However, the outflow is 
decelerated by the gravitational attraction of the galaxy 
hosting the AGN. Both effects increase with galaxy mass, 
resulting in competition. 



Fig. 2. — Final radius R of the outflows vs. galaxy mass Af ga i, 
for galaxies born in the redshift range 2bir = [3.1,3.0] (middle 
curve) and 2bir = [1.1,1.0] (bottom curve). Each dot in these 
curves represents one outflow. The top line shows a power-law 
approximation to these curves. 

These outflows do not reach equilibrium, and are still 
in the active- AGN phase or the post- AGN phase by the 
present time. The closer to the present the galaxies are 
born, the less time outflows have to propagate. 

Figure [H shows the final radius of the outflows vs. 
galaxy mass, for galaxies born in the redshift ranges 
Zbir = [3.1, 3.0] and Zhix = [1.1, 1.0]. More massive galax- 
ies produce outflows that travel larger distances. How- 



4.2. Total Masses in Simulation Volume 

Figure [3] shows the redshift evolution of the total bary- 
onic mass contained in galaxies hosting AGNs, the total 
mass of metals carried by outflows, and the enriched vol- 
ume fraction, inside the whole (182.6 Mpc) 3 simulation 
volume. The masses were computed by summing over 
AGNs with birth redshifts Zui > z. The top panel shows 
various total masses: the baryonic mass of all host galax- 
ies, ^ Mb, and the mass of metals transported by out- 
flows from all galaxies, ^2 -Mz.outj which includes metals 
carried from the central AGN, Mz,i, and metals car- 
ried from the non-central ISM, Mz,2 [see eq. Q]. The 
middle panel depicts the ratio of total metal mass over 
galaxy mass. Note that these total masses and their ra- 
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Fig. 4. — Total masses in the whole (182.6 Mpc) 3 simulation 
volume as a function of redshift, for AGNs with different bolomct- 
ric luminosity limits: the whole population, 10 8 < L/Lq < 10 14 
(solid), only the high-luminosity sources, 10 9 < L/Lq < 10 14 (dot- 
ted), only the low-luminosity sources, 10 8 < L/Lq < 10 9 (dashed). 
The top panel shows the total baryonic mass of the host galaxies 
HT) Mb) within the relevant luminosity range. The bottom panel 
shows the total mass of metals transported by outflows from those 
galaxies (J2 M z ,out)- 

tio do not depend on the outflow opening angle a. Thus 
the upper two panels are the same for runs A, B, and C. 
The bottom panel gives the fractional volume of the sim- 
ulation box enriched by the outflows in run C (a = 60°). 
This enriched volume fraction is the same as the solid 
curve in Figure 4 of Paper I, but is plotted here on a 
logarithmic scale. 

The total (dark matter + baryon) mass inside our 
simulation volume is Mbox = 2.32 x 10 17 M Q . Scaling 
by the cosmological baryon to total matter density ra- 
tio, this corresponds to a total baryonic mass in the 
box: M 6;box = M bax (Sl bfi /Sl ) = 3.84 x 1O 16 M . We 
can read out the values of masses from Figure |3J At 
the present epoch: the total baryonic mass in galax- 
ies, J2 M b( z = 0) = 3.18 x 10 14 M Q , and the total 
mass of metals available to be distributed to the IGM, 
E M z,out(2 = 0) = 1.18 x 1O 12 M . The ratio of total 
baryonic mass of galaxies over that of the simulation box 
is 0.00828. The total metal mass is a fraction 0.00370 of 
the baryonic mass in galaxies, and 3.07 x 10~ 5 of baryonic 
mass in the simulation box. Notice that only galaxies 
that host AGNs are included in our calculations. Galax- 
ies comprise of order 5% of the total mass in the uni- 
verse, hence galaxies hosting AGN represent a fraction 
0.00828/0.05 = 17% of all galaxies by mass. This is a di- 
rect consequence of the QLF we adopted for generating 
the AGN sources to populate our simulation volume and 
our assumption of the BH-galaxy mass relation ($2Tj). 



To study the relative contributions of low- and high- 
luminosity AGNs on metal enrichment, we divided the 
AGN population in two separate luminosity ranges, 
10 8 < L/Lq < 10 9 (low- luminosity or, faint AGNs) 
and 10 9 < L/Lq < 10 14 (high-luminosity or, bright 
AGNs) . Figure U shows the evolution of the baryonic 
mass in galaxies and the metal mass separately for low- 
and high-luminosity AGNs. We note that the total 
masses for the whole population (10 s < L/Lq < 10 14 ) 
is completely dominated by the high-luminosity sources. 
Hence the solid and dotted curves in Figure [4] almost 
coincide, and are not distinguishable in this plot. The 
low-luminosity AGNs (10 8 < L/Lq < 10 9 ) have a very 
small contribution to the total masses, at z — they 
comprise of 1.7% of the total galaxy baryonic mass and 
0.1% of the total metal mass. Summed over the entire 
population of fainter sources, the total mass of galax- 
ies is: 

Efaint( M gai) = 5 - 29 x 10 12 M Q , and the total 
mass of metals available to be distributed to the IGM: 
Efaint( M z,out) = 1-22 x 1O 9 M . The corresponding 
numbers for the brighter sources are: Ebright(-^gai) = 
3.12 x 10 14 M Q , Ebri g ht( M z,out) = 1.175 x 10 12 M Q . This 
point is relevant to explain the results about low- and 
high-luminosity AGNs we present in §4.31 and MA\ be- 
low. 

In Table [2] we list the contribution of AGNs in six dif- 
ferent bolometric luminosity ranges to the total baryonic 
mass of galaxies and total metal mass in the simulation 
volume. This is complementary to the information pre- 
sented in Figure 1 of Paper I. The first column indicates 
the luminosity range. Columns 2, 3, and 4 indicate the 
number of AGNs, the total baryonic mass of the host 
galaxies, and the total mass of metal, respectively, with 
the corresponding percentages in columns 5, 6, and 7. 
Even though low-luminosity AGNs dominate by num- 
bers, with 76.8% AGNs having L < 10 9 L Q and 94.0% 
AGNs having L < lO 10 ^©, it is the high-luminosity 
AGNs that contain most of the baryonic mass and pro- 
duce most of the metals. Only 1.6% of AGNs have a 
luminosity in the range 10 11 L Q < L < lO 13 ^©, but they 
contain 77.2% of the baryonic mass and produce 84.4% 
of the metals. 

4.3. Average Metallicity of the IGM 

Using equations (0 and ([6]), we calculated the metal 
content and metallicity caused by the AGN outflows, in 
a Nfi — 256 3 grid, for each run. Figure [5] shows the 
maximum metallicity [0/H] cc n iinax in the computational 
volume versus redshift, for outflows with various opening 
angles. The maximum metallicity increases with time up 
to redshift z = 2. As new AGNs are formed and produce 
outflows, these outflows overlap with the ones produced 
by earlier AGNs, causing an increase in metallicity. After 
z = 2, the value [0/H] ce n imax levels-off around 0, and 
the variations with z are no longer monotonic. At that 
point, many outflows have reached the post- AGN phase. 
They are still expanding, but no additional metals are 
being added, so the metals contained in the outflows are 
diluted. This effect competes with the addition of metals 
by late-forming AGNs that are still in the active-AGN 
phase. 

We compute the mass-weighted and volume-weighted 
average metallicities induced in the IGM in our Simula- 
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TABLE 2 

Galaxy and metal mass for AGNs with different bolometric luminosities 



] NAGN (Y, M b)Wo\ (E(Af Z ,out) [Mq] %N %M h %M Z ,out 



10 8 


- 10 9 


1178571 


5.29 x 10 12 


1.22 x 10 9 


76.8 


1.7 


0.1 


10 9 


- 10 10 


264676 


1.22 x 10 13 


4.86 x 10 9 


17.2 


3.9 


0.4 


10 10 


-10 11 


68655 


3.24 x 10 13 


3.45 x 10 10 


4.5 


10.2 


2.9 


10 11 


-10 12 


19372 


9.10 x 10 13 


2.49 x 10 11 


1.3 


28.7 


21.2 


10 12 


-10 13 


3978 


1.54 x 10 14 


7.44 x 10 11 


0.3 


48.5 


63.2 


10 13 


-10 14 


110 


2.26 x 10 13 


1.43 x 10 11 


0.0 


7.1 


12.2 
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Fig. 5. — Maximum metallicity attained by a cell in the whole 
simulation volume, as a function of redshift, for different outflow 
opening angles: a = 60° (solid), 120° (dotted), 180° (dashed). 

tions. The averages are computed by considering only 
the volumes which are enriched by the outflows. Hence, 
only cells which contain enriched particles (among Ng 
total in the whole simulation box, §2.4|) are included. In 
computing the metallicity [O/H], the log 10 is taken after 
doing the averages over the relevant volumes. Mathemat- 
ically, the global volume-weighted average metallicity in 
the enriched volumes is expressed as: 

[O/H] 

volAvg — 

[(O/H)] 

volume 

j2 ( ^) I ( ^ 

died colls ^ " H ' ccl1 ' ^ ™ H 



log 



10 



enriche 



Number of enriched cells 



(7) 



(Note: this is effectively a volume-weighted average, 
since all cells have the same volume). The global mass- 
weighted average metallicity in the enriched volumes is 
given by: 



[0/H] raassA vg = [(O/H)], 



logi 



enriched 



lied colls H / cell / V H / ( 



J2 /'cell 
enriched cells 



(8) 



For the solar oxyge n abundance we use: 12 + 

logio( n o/ n H)o = 8 - 93 (|Anders fc Grevessdll989T ). 
The redshift evolution of [0/H] vo i Avg and [0/H] massAvg 

in the enriched volumes of the simulation box are pre- 
sented in Figure [B] for AGN outflows with various open- 
ing angles, and Figure [7] for AGN populations with dif- 
ferent bolometric luminosity limits. We note that the 
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Fig. 6. — Average metallicity (volume- weighted average in the 
top panel, and mass-weighted average in the bottom panel) in the 
enriched volumes of the simulation box, as a function of redshift, 
for different outflow opening angles: a = 60° (solid), 120° (dotted), 
180° (dashed). 

values of volume-weighted and mass-weighted averages 
in a simulation are quite similar. In runs A, B, C, and 
Cl > [0/ H ]voiAv g is sli S htl y la rg er than [0/H] massAvg at 
high-z, while after a certain redshift (z ~ 1.1 in run A, 
0.7 in B, 0.3 in C and CI) [0/H] massAvg becomes higher. 
A reverse trend is seen in run C2, where [0/H] massAvg 
is greater at high-z, and [0/H] volAvg takes over after 
z ~ 0.5. 

The whole population of AGNs (having bolometric lu- 
minosities between 10 s < L/L & < 10 14 ) with most- 
anisotropic outflows (opening angle a — 60°, run C) in- 
duces a volume-averaged IGM metallicity of [O/H] f=s 
—4.9 at z = 5.5, peaking to [O/H] w —2.6 between 
z ~ 1 — 1.5, and then leveling off (slightly decreasing) to 
[O/H] ss -2.9 at the present epoch. Run B (a = 120°) 
and the isotropic outflow case (a = 180°, run A) produce 
quite similar average results: starting with [O/H] « —5.1 
at z = 5.5, then rising gradually, and lying more or 
less flat at [O/H] w -2.8 between z = 2 - 0. There 
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Fig. 7. — Average metallicity (volume-weighted average in the 
top panel, and mass-weighted average in the bottom panel) in the 
enriched volumes of the simulation box, as a function of redshift, 
for different bolometric luminosity limits: the whole population, 

10 8 < L/Lq < 10 14 (solid), only the high-luminosity sources, 

10 9 < L/Lq < 10 14 (dotted), only the low-luminosity sources, 
10 s < L/Lq < 10 9 (dashed). 

is clearly a change of slope around z — 2, where the 
mass-averaged metallicity and volume-average metallic- 
ity level off. Again, this can be attributed to the fact 
that many outflows have reached the post-AGN phase. 
The resulting spreading of the metals compensates for 
the injection of new metals by late-forming AGNs. 

The variations of these results between runs A, B, and 
C are small. Outflows with different opening angles pro- 
duce comparable average metallicities in the IGM. How- 
ever, we find that more anisotropic outflows (a = 60°) 
causes higher average metallicity at high-z, and the 
trend reverses at low-z, with more isotropic outflows 
(a = 180°) giving larger metallicity values. At z = 0, 
the average metallicities are essentially the same for all 
opening angles. For the volume- aver age metallicity, this 
was expected, since the amount of metals deposited by 
AGNs is the same for all opening angles, and the en- 
riched volume fraction is also the same (Paper I). For 
the mass-averaged metallicity, the extent of overlaps be- 
tween outflows could have made a difference, but this is 
clearly not the case. 

The high-luminosity AGNs (10 9 < L/Lq < 10 14 , 
run CI) produce slightly higher metallicities at high-z; 
starting with [O/H] ks —4.6 at z = 5.5, and converge 
to similar values as the whole population (run C) by 
z ~ 3. The contribution of the low-luminosity AGNs 
(10 8 < L/Lq < 10 9 ) to the resultant IGM metallicity 
is much lower, 2 — 3 orders orders of magnitude smaller 
than the luminous sources. Run C2 produces an aver- 
age metallicity of [O/H] w —7 at z = 5.5, which only 



increases to [O/H] « —6.4 at the present epoch. We 
therefore conclude that the metallicity induced in the 
IGM by anisotropic AGN outflows is greatly dominated 
by sources having bolometric luminosity L > 10 9 Lq, 
sources with 10 8 < L/Lq < 10 9 have a negligible contri- 
bution. 

We compare our results with other observational and 
theoretical studies. As detailed in SfTJ C IV and O VI ob- 
servations show that the average met allicity of the IGM 
is Z > 10~ 3 - 10~ 2 ^ at z ~ 2 - 3 dCowie et al.lfl995l: 
Songaila fc Cowid 119961; iSongailal 119971 : iCarswell et al.l 
20021: iSimcoe et alll2004l ). with few indications that an 
IGM me tallicity of ~ 10~ 4 Zq is already in pl ace at 
z = 5 (jSongailal [200TI : IRvan- Weber et al.l 1200!):. At 
low redshifts, a mea n cosmic metallicity of Z > 10~ 2 — 
10~\Zn is detected (iHurles fc Tvtlerl fl996t iTripp et al.l 
12001 IDanforth fc Shullll207)5n ~ 

Several studies and cosmological numerical simula- 
tions have shown that ejection of metals via winds 
from galaxi es can reprodu c e the observed IGM metal 
properties ( Tcgm ark et al.l 119931 : lAguirre et al.l 120011 : 
iScannapieco et al.l 120061 ) .These simulations have shown 
that outflows driven by SNe or starbursts can produced 
IGM mean metallicity values consistent with observa- 
tion. In these simulations, most of the outflows orig- 
inate from dwarf and low-mass gala xies, and start at 
an ea rly epoch (6 < z < 15). IScannapieco et all 
(2002) found that most of the enrichment occurred rel- 
atively early, with mass-average d cosmological metal- 
licity values 10~ 3 — 10~ 15 Z Q . lOppenheimer fc Pavel 
(2006) obtained a global average metal mass fraction of 
log(Z/Z ) ~ —2 at z ~ 2. We found similar values for 
the volume-weighted and m ass-weighted average, while 
Oppenhcimer & Dave (2006) obtained a mass-weighted 
average larger than the volume- weigh ted average. This is 
because lOppenheimer fc Pavel (|2006l ) averaged over the 
entire computational box, whereas we average only over 
the IGM volumes which are enriched in metals. 

The values of IGM metallicities we computed are some- 
what lower compared to these observations and simula- 
tions. Our simulations produced an average IGM metal- 
licity of [O/H] ~ — 5 at z = 5.5, then rising gradually, and 
lying relatively flat at [O/H] —2.8 between z — 2 — 0. 
Comparing these results with observations, we estimate 
that AGN outflows contribute more than 10% of the IGM 
metallicity at z = 5, between 16% and 100% at z = 2, 
and between 2% and 16% at z = 0. The biggest factor 
coming to play here is the amount of metals in the galax- 
ies which is available to be ejected to the IGM by the 
outflows. The total galaxy baryonic mass of AGN hosts 
(galaxy population obtained from the QLF) is 0.0083 of 
the total baryons in the simulation box at z = 0. AGNs 
are only hosted in a fraction of galaxies of the Universe, 
and our model tracks only the metals generated in the 
AGN host galaxies. Since we do not take into account 
metals generated in non-AGN host galaxies, over the cos- 
mic epochs, the resulting IGM metallicities we obtained 
are smaller than (~ 10% of) the observed values. 

Some other simulation-based studies indicate that less- 
massive galaxies dominate in enriching the IGM at high- 
z, and ejection by AGN outflows accounts for only a 
fraction of the IGM metallicity. iCen fc Brvanl ()200lD 
stressed that dwarf and subdwarf galaxies with masses 
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10 6 5 — 1O 9 M are largely responsible for the metal pollu- 
tion of the I GM (~ 10~ 2 Z p) , as o bserved in Ly-a clouds) 
at z ~ 3. Thack er" et al.l {2002) showed that outflows 
from galaxies with total halo masses < 1O 1O M enrich 
~ 20% of the simulation volume with a mean metallicity 
of 0.3% solar at z = 4. 

Few direct observations of powerful outflows from AGN 
have been used to constrain the cosmological contribu- 
tion of AGN outflows to enriching the IGM. Extend- 
ing the observed properties of the rad io galaxy MRC 
1138-2 62 using cosmological constraints, Nesvadba et al.l 
(2006) estimated that powerful AGN feedback ejects 
~ (1 - 30) x 1O 4 Z/Z M Mpc- 3 of metals. So AGN 
feedback from massive galaxies may account for only few 
to ~ 20% of all metals in the IGM at z ~ 2. Com- 
paring the relative impact o f AGN - vs. starburst-driven 
winds. iNesvadba h Lehnertl (|2007D did some preliminary 
calculations using the comoving number density of high- 
z quasars, corrected for their duty cycle and energy in- 
jected. They found that the metal mass expelled by 
AGN winds into the IGM is smaller than that by star- 
burst winds, concluding that AGN winds are likely not 
the dominant mode of IGM metal enrichment through 
cosmic time, however they could be a significant contrib- 
utor. Our results are consistent with such works that 
propagation of metals by AGN outflows do not account 
for 100% of the observed metals in the IGM. We con- 
clude that ejection of metals from (those generated in) 
AGN host galaxies by AGN outflows enriches the IGM to 
> 10 — 20% of the observed values, the number dependent 
on redshift. 

This could give the impression that metal enrich- 
ment by AGN-driven outflows is a minor effect com- 
pared to metal enrichment by SNe-driven outflows. How- 
ever, one must take into account the enriched volume 
fraction, which is of order 10 — 20% for SNe-driven 
outflows dScannapieco k. Broadhurst 2001; Madau et al.l 



120011: iBertone et all 120051: iPieri et all 120071). and over 



80% for AGN-driven outflows (|Levine fc Gnedml 120051: 
Paper I). The final radii of AGN outflows shown here 
in Figures [1] and [5] are s ignificantly large r than the ones 
for SNe-driven outflows (jPieri et al.ll2007l Figure 8). The 
SNe-driven outflows do not travel far from the cosmologi- 
cal structures from where they originate, and they are not 
powerful enough to reach the deepest voids. Thus, AGNs 
are likely one of the most important contributors to metal 
enrichment of the IGM, especially the low-density regions 
because of their ability to trigger powerful outflows. 

4.4. Enriched Volume above Metallicity Limit 

We calculated the volume fractions of the simula- 
tion box which are enriched above metallicity limits of 
[O/H] = —4, —3, and —2.5. In practice, we count the 
number of cells in the computational volume which have 
[0/H] coll > -4, [0/H] cell > -3, and [0/H] coll > -2.5 
and divide by the total number (iVf) of cells. Figure [8] 
shows the redshift evolution of the fractional volumes 
above these metallicity limits for AGN outflows with var- 
ious opening angles. At high-z, more anisotropic out- 
flows (lower values of a) enrich slightly larger volume 
fractions. The trend reverses at low-z, where run B 
(a = 120°) and the isotropic outflow case (a = 180°, 
run A) enrich similar volume fractions, somewhat higher 




Fig. 8. — Volume fractions enriched above metallicity limits: 
[O/H] > —2.5 in the top panel, [O/H] > —3 in the middle panel, 
[O/H] > —4 in the bottom panel, as a function of redshift, for 
different outflow opening angles: a = 60° (solid), 120° (dotted), 
180° (dashed). 
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Fig. 9. — Volume fractions enriched above a certain metallic- 
ity limit, ^/([O/H] > [O/HJjj it ), at the present epoch, for differ- 
ent outflow opening angles: a = 60° (solid), 120° (dotted), 180° 
(dashed). 

than the most- anisotropic outflows (a = 60°, run C). 
These enriched volume fractions are small at z > 3, then 
rises rapidly. We found that the volumes enriched above 
these metallicity limits by the high-luminosity sources 
(10 9 < L/Lq < 10 14 , run CI) are < 1% lower than that 
by the whole population of AGNs (10 8 < L/L Q < 10 14 , 
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TABLE 3 

Average Metallicity and Enriched Volume Fraction atz = 



Run a (°) [Q/H] volAvg y frac (> -4) V bac (> -3) V frac (> -2.5) 



A 


180 


-2.830 


0.453 


0.234 


0.098 


13 


120 


-2.877 


0.449 


0.225 


0.087 


C 


60 


-2.960 


0.345 


0.146 


0.056 


CI 


60 


-2.934 


0.344 


0.146 


0.056 


C2 


60 


-6.424 


0.000 


0.000 


0.000 



run C). This is not distinguishable in this plotting scale, 
so we do not show them separately. 

In Table[3]we list some of the main results presented in 
this section and the previous one. The first 2 columns are 
the same as in TableQ] The third column gives the result- 
ing volume- weighted average metallicity ([0/H] volAvcr ) at 
the present epoch. The fourth, fifth, and sixth columns 
give the volume fraction enriched above metallicity of 
[O/H] = —4, —3, and —2.5 at the present epoch, re- 
spectively. While the average metallicity is quite insensi- 
tive to the opening angle (varying from —2.83 to —2.96), 
There is a significant dependence of the volume fraction 
on the opening angle, with larger opening angles result- 
ing in larger values. At larger opening angles, there is 
more overlap between outflows, resulting in larger val- 
ues of the metallicity in the regions where overlap oc- 
curs. The low-luminosity AGNs (10 s < L/L Q < 10 9 , 
run C2) do not enrich any volume to any of these metal- 
licity limits ([O/H] > —4) at any time. This is not sur- 
prising, since the average metallicity it produces is only 
[O/H] w -6.4 at present. 

We compa re these results with some related contempo- 
rary studies. iScannapieco et all (|2002l ) performed Monte 
Carlo cosmological simulations to track the time evo- 
lution of SN-driven metal-enriched outflows from early 
galaxies, finding that up to 30% volume of the IGM is en- 
riched to above 10 _3 Zq at z — 3, making the enrichment 
biased t o the areas near the st a rburst ing galaxies them- 
selves. lOppenheimer fc Dave! (|2006t ) found that their 
resulting volume fraction having [Z/H] > — 3 reaches 
~ 10% at z = 2. Qp penheimer et all (|2009( ) calculated 
that the volume filling factor of metals increases between 
z = 8 and z = 5, reaching ~ 1% for Z > 10~ 3 Z Q by 
z = 5. The volumes we obtained enriched above given 
metallicity limits are smaller than these studies, for rea- 
sons same as that discussed in §4.31 We only eject the 
metals produced in AGN host galaxies. At high redshifts, 
AGN-driven outflows enrich smaller volume fraction of 
the IGM above a certain metallicity than SNe-driven out- 
flows which transport metals from less-massive galaxies 
as well. 

The continuous behavior of the enriched IGM volumes 
as a function of metallicity is presented in Figure [SJ 
It shows the fractional volume enriched above a cer- 
tain metallicity limit at z = 0, caused by AGN-driven 
outflows with various opening angles. We see that the 
more isotropic outflows (a = 120° and 180°) enrich al- 
most the same volume fractions above a given metallic- 
ity, 5 — 10% higher than the most-anisotropic outflows 
(a = 60°). More than 60% of the volume is enriched 
above [O/H] > —6, but only and 1 — 2% of the volume 
is enriched above [O/H] > —2. The curves we obtained 
in Figure |H] are qualitatively similar to the one shown in 



Figure 5 of lThacker et al.l (|2002l ) . Even though the scales 
and redshift are different (and they consider SNe-driven 
outflows), the trends are similar. At high abundances, 
there is a rapid increase in enriched volume fraction with 
decreasing metallicities, while at lower abundances the 
enriched volume fraction tends to get flatter. 

4.5. Density - Enrichment Correlation 

We compute the mean metallicity of the IGM at a given 
density, [0/H] p , within the range —3 < \og{p/p) < 4, 
where p — iVl^H^(l + z) 3 /8irG is the mean density of 
our simulation box at redshift z. We divide the density 
range in bins of size 0.1 dex, and identify the cells which 
are enriched and whose total matter density p ce \\ occur 
within each bin (see q2.4j) . We average the gas density 
Pg = Pccii^f>,o/^o and metal density Pmetai within each 
bin to get the bin-averaged values p g ,bin, Pmctahbin, and 
■Zcdi.bin = Pmotai,bin/p g ,bin- We then calculate the mean 
metallicity [O/H],, = log(n /n H ) (9 - log(n o /n H ) inside 
each bin, where log(7i /n H ) p is calculated using equa- 
tion ©. 

The evolution of the resulting metallicity as a function 
of IGM density (total density, including both baryons 
and dark matter) is shown in Figure [101 at redshifts 
z = 4, 2, and 0, in the columns from left to right. 
As the simulation proceeds, overdense regions collapse 
while voids become deeper, increasing the range of den- 
sities present in the simulation volume. Also, with time 
the population of outflows cover a larger fraction of the 
simulation volume, causing densities of a larger range 
(denser filaments and less dense void regions) to become 
enriched. As a result, the range of densities (cc-axis) cov- 
ered by each curve increases from z = 4 to z — 0. At 
a given total density, the metal density (p me tal,bin) top 
panels) increases with decreasing redshift. This happens 
because with time new outflows are born, which continu- 
ously add more metals to the IGM volumes. The middle 
panels show the gas density (p g) bin) of the enriched IGM 
volumes, which by construction follow the total density 
at all times (we assume p g = f^op/ilo)- 

The bottom panels of Figure [TO] show the metallic- 
ity [0/H] p , which grows with time steadily (at a given 
IGM density). This is consis tent with the trend found by 
lOppenheimer fc Pavel (|2006h in their Figure 10. We find 
some interesting trends in the metallicity-density corre- 
lations as well as some differences between outflows with 
different opening angles. At higher redshifts (z = 4 and 
2) the underdense regions [log(p/p) < 0] are enriched 
to higher metallicities, and the induced metallicity de- 
creases with increasing IGM density. This trend is more 
prominent with increasing anisotropy of the outflows. It 
is very clear in the bottom-left and bottom-middle panels 
in the plot for run C (a = 60°). At z = 2, the metallicity 
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Fig. 10. — Metallicity induced by AGN outflows as a function of the total density of the IGM, at redshifts 2 = 4 (left column), z = 2 
(middle column) and 2 = (right column), top: metal density, middle: gas density, bottom: metallicity (logarithm of ratio of the top to the 
middle panel, minus logarithm of the ratio for the solar values). The line types indicate different outflow opening angles: a = 60° (solid), 
120° (dotted), 180° (dashed). 

still decreases with increasing density in overdense re- 
gions [log(p/p) > 0] for all opening angles. At z = 0, the 
metallicity is essentially flat over all IGM densities, for 
outflows with any opening angle. At that redshift, the 
enriched volume fraction exceeds 80% (Paper I), hence 
all regions, at all densities, are enriched in metals, except 
the deepest voids [\og(p/p) < —2.4] located far from any 
AGN. Initially, more anisotropic outflows preferentially 
enrich low-density regions, but this trend gets eventually 
washed-out as the enriched volume fraction approaches 
unity. 

Our trend of larger metallicity induced in the IGM at 
low densities for the more anisot ropic outflows at high- 
z is o pposite from the results of IQppenheimer fc Pavel 
(2006), who found a steady increase of metallicity with 
overdensity. The reason for this discrepancy is the 
anisotropy of the outflows in our model. The AGN out- 
flows propagate along the direction of least resistance 
around a density peak ( fl2.1[) . and preferentially enrich 
the low-density regions to higher metallicities at early 
epochs. In our simulations, the isotropic outflow case 
(a = 180°) produces uniform metallicity at all densities 
being enriched. 

As we have detailed in [JT] metal enrichment of the IGM 
is found to be highly inhomogeneous, strongly dependent 
on local density and redshift. The first detection of O VI 
in the low-dens i ty IG M at high redshift was reported 
by iSchave et al.l ((20001 ) , who performed a pixel- by-pixel 
search for O VI absorption in quasar spectra over the red- 



shift range z = 2 - 4.5. They detected O VI at 2 < z < 3 
down to r HI ~ 10 _1 , showing that the IGM is enriched 
down to much lower overdensities, and did not detect 
O VI at z > 3. Enrichment by powerful AGN outflows at 
high-z (following our model prescriptions) can naturally 
explain observations of such metal-enriched underdense 
regions. Providing s upport to such a AGN- driven IGM 
enrichment scenario, iKhalatvan et al.l ((2008) found, us- 
ing SPH simulations, that without AGN feedback metals 
are confined to vicinities of galaxies, underestimating the 
observed metallicity of the IGM at overdensities < 10 
(the underdense IGM). 



4.6. Diffuse and Dense IGM 

To complete the analysis of our simulations, we divide 
the whole simulation volume into two components: the 
diffuse IGM having density p/p < 100 and the dense 
iGM with density p/p > 100. We compute the mass- 
weighted average metallicity in the enriched volumes for 
the diffuse and dense IGM using, 
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Fig. 11. — Mass fractions and metallicity for the diffuse (p/p < 
100, upper curves) and the dense (p/p > 100, lower curves) regions 
of the simulation volume, as a function of redshift. Total mass 
fraction in the top panel, metal mass fraction in the middle panel, 
and mass- weighted average metallicity within the enriched volumes 
in the bottom panel). The line types correspond to different outflow 
opening angles: a = 60° (solid), 120° (dotted), 180° (dashed), as 
indicated. 
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Figure [TT] shows the redshift evolution of various quan- 
tities (mass fraction, metal fraction, and mass-weighted 
average metallicity within the enriched volumes) in the 
diffuse and dense IGM, for different outflow opening an- 
gles. In each panel, the upper set of curves is for the 
diffuse volumes (p/p < 100), and the lower set of curves 
denote the dense volumes (p/p > 100). 

As time goes on, the mass fraction in the dense vol- 
umes increases as condensed structures grow within the 
simulation volume, by accreting matter from the diffuse 
IGM, whose mass fraction then decreases. We find a 
similar trend of the metal fraction, the diffuse regions 
containing more metals than the dense regions. Con- 
sequently, as seen in the lower panel of Figure [TTJ the 
average metallicity of the diffuse volumes is larger than 
that of the dense volumes. These trends are contrary 



to the ones found bv lQppenheimer fc Pavel ([2006D (their 
Fig. 7). There are two reasons for this discrepancy. First, 
as we explained before, they have included all galaxies in 
their simulations, while we only include AGN-host galax- 
ies. Second, they consider SNe-driven outflows only while 
we consider AGN-driven outflows. SNe-driven outflows 
achieve enriched volume fractions of order 20% compared 
to over 80% for AGN-driven outflows, which are therefore 
much more efficient in enriching low-density regions. 

4.7. Different Initial Location Scheme: More Luminous 
AGN in Denser Peaks 

The results presented in the previous subsections are 
based on the 5 simulation runs done using the same fidu- 
cial AGN initial location prescription as in our Paper I, 
as mentioned in ' i'Z.W Restating the method: we spatially 
locate the new AGNs born during a timestep interval at 
the local density peaks, selected randomly among all the 
peaks in the cosmological volume at that timestep of the 
simulation. (Runs H and I of our Paper I considered 
density bias and clustering in locating the AGNs.) 

In order to test the robustness of this fiducial AGN ini- 
tial location scheme w.r.t. AGN-luminosity-peak-dcnsity 
correlations, we performed a test simulation run by im- 
plementing a new prescription. In this run, after se- 
lecting the peaks at each timestep to locate the new 
AGNs, the peaks were sorted according to their densities 
and the AGNs were sorted according to their luminosi- 
ties. Then each AGN was assigned a peak such that the 
highest-luminosity AGN is located in the highest-density 
peak, and so on. This method was repeated at all the 
timesteps. 

The results of this test run came out to be almost 
similar to our corresponding run with the fiducial loca- 
tion scheme. Figure shows a comparison of the re- 
sults between the two runs. The top-left panel shows 
that the mass-weighted average metallicity is undistin- 
guishable, and the new run produced a slightly higher 
volume-weighted average metallicity. The volume frac- 
tions enriched above metallicity limits are undistinguish- 
able (from the top-right and bottom- right panels). The 
quantity that showed the relatively largest difference is 
the metal mass fraction in the diffuse and the dense re- 
gions of the simulation volume, shown in the bottom-left 
panel. At the same time, the extent of this slight differ- 
ence reduces at lower redshifts (z < 2). 

Overall there is a ~ 1% difference. This shows that 
our results are not sensitive to correlating the AGN- 
luminosity with the peak-density in our AGN initial lo- 
cation prescription. In our model, at every timestep, we 
use the observed QLF to calculate the number of AGNs 
born during that timestep, and we locate these AGNs in 
density peaks that are available at that timestep. Such 
a technical implementation averages out any effect aris- 
ing from correlating the AGN-luminosity with the peak- 
density. 

5. SUMMARY AND CONCLUSION 

We have computed the metallicity of the IGM induced 
by outflows of the cosmological population of AGNs over 
the age of the Universe. In Paper I, we developed a semi- 
analytical model for anisotropic AGN outflows expand- 
ing along the path of least resistance, and implemented it 
into ./V-body cosmological simulations. In this paper, we 
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Fig. 12. — Comparison of results with the new AGN initial location prescription vs. the fiducial scheme. The older results (run C) are 
shown as solid lines. The new run is shown by dotted lines, lying very close to the old run in each case. Top-left: Volume-weighted average 
(upper curves), and mass- weighted average (lower curves) metallicity, in the enriched volumes of the simulation box. Top-right: Curves, 
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Bottom- right: Volume fraction enriched above a certain metallicity limit at the present epoch. 



have incorporated a detailed metal enrichment technique 
in our AGN outflow model from Paper I, and have simu- 
lated the metal-enrichment of the IGM. Our enrichment 
prescription consists of estimating the masses of metals 
transported with the outflows from the host galaxies, and 
distributing the metals to the IGM volumes being inter- 
cepted by the outflows. 

During their expansion, the AGN outflows carry with 
them metals generated within the host galaxy (by stars) , 
and spread the metals to the broader IGM. We distin- 
guish 2 regions of stellar populations as the sources of 
metal generation. For the stars near the AGN, we con- 
sider a central metal abundance of 5Zq (consistent with 
observed super-Solar metallicity of broad-line regions). 
We further take into account that AGN outflows can en- 
train some enriched ISM of the greater galaxy and spread 
it to the IGM. For stars located away from the central 
AGN within the rest of the galaxy, we use the redshift- 
dependent mass-metallicity relation of galaxies derived 
from observations. Enrichment of the IGM is simulated 
by distributing the metals transported by the outflows to 
the intergalactic gas intercepted by them. Each outflow 
is assumed to impart all of the stripped metals from its 
host galaxy to the IGM over the active-AGN lifetime, 
depositing metals at an uniform rate in time. 

Using this algorithm, we simulated the propagation 
of AGN-driven outflows carrying metals, along with the 



growth of large-scale structures, in a cosmological vol- 
ume of size (128 h~ x Mpc ) 3 in a ACDM universe, and 
analyzed the resulting metal-enrichment of the IGM. We 
performed a series of 5 simulations by varying the out- 
flow opening angle and the AGN population limiting lu- 
minosities. The main results and conclusions from our 
study arc in the following. 

1) In our simulation box of comoving volume (128 h~ x 
Mpc) 3 , the total baryonic mass of all the AGN host 
galaxies is 3.18 x 10 14 M Q , and the total mass of met- 
als available to be distributed to the IGM (by the out- 
flows) is 1.18 x 10 12 M Q , at the present epoch. The ra- 
tio of total baryonic mass of galaxies over that of the 
simulation box is 0.00828. The total metal mass is a 
fraction 0.00370 of the baryonic mass in galaxies, and 
3.07 x 10~ 5 of the baryonic mass simulation box. The 
masses of the whole AGN population (10 8 < L/L & < 
10 14 ) is completely dominated by the high-luminosity 
sources (10 9 < L/L Q < 10 14 ). The low-luminosity AGNs 
(10 8 < L/Lq < 10 9 ) have a very small contribution to 
the total masses, at z = they comprise of 1.7% of the to- 
tal galaxy mass and 0.1% of the total metal mass. Even 
though high-luminosity AGNs are much less numerous 
than low-luminosity ones, they really dominate the mass 
budget, and form the dominant source of metals to be 
distributed to the IGM. 
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2) The values of volume-weighted and mass-weighted 
average metallicity induced in the IGM are quite similar, 
in each simulation. The whole population of AGNs (hav- 
ing bolometric luminosities between 10 8 < L/Lq < 10 14 ) 
induces an average IGM metallicity of [O/H] ss —5 at 
z = 5.5, which then rises gradually, and flattens at 
a value [O/H] ss —2.8 between z = 2 — 0. Outflows 
with different opening angles produce comparable aver- 
age metallicities in the IGM. However, more anisotropic 
outflows (a = 60°) causes slightly higher average metal- 
licity at high-z, and the trend reverses at low-z, with 
more isotropic outflows (a = 180°) giving larger metal- 
licity values. The high-luminosity AGNs (10 9 < L/L Q < 
10 14 ) produce slightly higher metallicities at high red- 
shift, starting with [O/H] w —4.6 at z = 5.5, and con- 
verging to similar values as the whole population by 
z ~ 3. The contribution of the low- luminosity AGNs 
(10 s < L/L Q < 10 9 ) to the resultant IGM metallicity 
is much lower, 2 — 3 orders of magnitude smaller than 
the luminous sources. We conclude that the metallic- 
ity induced in the IGM by anisotropic AGN outflows is 
greatly dominated by sources having bolometric luminos- 
ity L > 10 9 L Q , sources with 10 8 < L/L Q < 10 9 having 
a negligible contribution. 

3) Observational studies indicate that the average 
metallicity of th e IGM is Z > 10~ 3 - lO" 2 ^ at 
z ~ 2 — 3 (e.g., iSongaila fc Cowiel 119961 : iSimcoe et al.l 
2004), with few indications that an IGM metallicity 
of ~ 10~ 4 ZfD is already in p lace at z — 5 (jSongailal 
2001; Rv an- Weber et alj|2009ft . At low redshifts, a mean 
cosmic metallicity of Z > 10~ 2 - \0~ 1 Zp ) is detected 
([Tripp et all [20021 : IPanforth fc Shul [20051) . Some cos- 
molq g ical numerical simulations (e . g . , iScannapieco et all 
120021 : lOppenheimer fc Pavel &im ) have deduced IGM 
mean metallicity values consistent with observations, 
taking into account metals ejected by outflows driven 
by SNe or starbursts, most of which come from low- 
mass galaxies, starting from early epochs (6 < z < 
15). Our computed results of average IGM metallicity 
([O/H] = -5 at z = 5.5, then rising to [O/H] = -2.8 
at z = 2 — 0) are somewhat smaller compared to these 
observations and simulations. The reason is that our 
model tracks only the metals produced and ejected from 
galaxies hosting AGNs, which is a fraction of the cosmo- 
logical population of galaxies in the Universe, (we do not 
take into account metals generated in non-AGN galax- 
ies), and that cannot account for 100% of the metals 
observed in the IGM. We conclude that ejection of met- 
als from AGN host galaxies by AGN outflows enriches 
the IGM to > 10 - 20% of the observed values, the 
number dependent on rcdshift. Some other simulation- 
based work indicates that less-massi ve galaxies dominate 
the enrichment the IGM at high-z (|Cen fc Brvanl 120011 : 
iThacker et al.ll2~002t ). supporting the scenario that ejec- 
tion by AGN outflows accounts for only a fraction of 
the IGM metallicity. Fro m direct observat i ons of pow- 
erful outflows from AGN, iNesvadba et all (|2006fl esti- 
mated that powerful AGN feedback from massive galax- 
ies may account for only few to ~ 20% of all metals in the 
IGM at z ~ 2. Comparing the relative impact of AGN - 
vs. starburst-driven winds, Ncsvadba & Lchncrt (2007) 
concluded that AGN winds are likely not the dominant 
mode of IGM metal enrichment through cosmic time, 



however they could be a significant contributor. 

4) At high-z, more anisotropic outflows (lower values 
of a) enrich slightly larger volume fractions. The trend 
reverses at low-z, where the runs with a = 120° and 
180° enrich similar volume fractions, somewhat higher 
than the most-anisotropic outflows (a = 60°). These 
enriched IGM volume fractions above certain metallic- 
ity limits are small at z > 3, then rises rapidly to the 
following values at the present epoch: 5.8 — 10% of the 
volume enriched to [O/H] > -2.5, 14.5 - 24% volume to 
[O/H] > -3, and 33.8-45% volume to [O/H] > -4. The 
more isotropic outflows (a = 120° and 180°) are seen to 
enrich almost the same volume fractions above a given 
metallicity, 5 — 10% higher than the most-anisotropic 
outflows (a = 60°). At z = 0, > 60% of the volume 
is enriched to [O/H] > —6, and 1 — 2% volume has 
[O/H] > —2. Our results of fractional volumes enriched 
as a function of metallic ity limits are consistent with that 
of IThacker et all (|2002T l. 

5) Examining the correlations of the induced enrich- 
ment with the IGM density, we find that the metallic- 
ity [0/H] p grows with time steadily (at a given IGM 
density). This is consisten t with the trend found by 
lOppenheimer fc Pavel (|2006h . At higher redshifts (z > 
2), for more anisotropic outflows, there is a prominent 
gradient of the induced enrichment, the metallicity de- 
creasing with increasing IGM density, making the un- 
derdense (p < p) regions enriched to higher metallicities. 
This trend is more prominent with increasing anisotropy 
of the outflows. The metallicity gradient (vs. density) 
decreases as the outflows become more isotropic (larger 
values of a). At z = 0, the metallicity is more or less flat 
over all IGM densities, for outflows with any opening 
angle. 

In Paper I, we had found that increasingly anisotropic 
AGN outflows preferentially enrich lower-density vol- 
umes of the IGM. In the present work, we complement 
these results by showing that more anisotropic outflows 
enrich the underdense IGM to higher metallicities at 
high-z. Evidence of substantial m etallicity in underde nse 
regions of the IGM at z - 4 (e.g- lSchaye et al.ll2003h re- 
quires a strong mechanism of spreading metals widely 
into low-density regions, at early cosmic epochs. SNe- 
driven outflows reach enriched volume fractions of order 
20% only, and therefore cannot enrich low-density re- 
gions easily. Enrichment by powerful AGN outflows at 
high-z (following our model prescriptions) can naturally 
explain observations of such metal-enriched underdense 
regions. 

6) We found that that our results are not sensitive to 
correlating the AGN-luminosity with the peak-density in 
our AGN initial location prescription. In other words, in 
our model when locating the AGNs born at a timestep in 
the density peaks available at that timestep, it does not 
make a diflcrcncc if the more luminous AGNs are located 
in the denser peaks, or whether the AGNs are located in 
the available peaks randomly. 

7) We found that the anisotropy of the AGN outflows 
affect most of the final results only marginally, except 
at highcr-z. Outflows with different opening angles pro- 
duced comparable volume-weighted and mass-weighted 
average metallicities in the IGM at all redshifts. The 
IGM volume fractions enriched above metallicity limits 
is larger for more anisotropic outflows at high-z, while 
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at z < 1.5 more isotropic outflows enriched a larger frac- routine that calculates the direction of least resistance 

tion. At z > 2, more anisotropic outflows enriched the was written by Cedric Grenon. All calculations were per- 

underdense IGM to higher metallicities. formed at the Laboratoire d'astrophysique numerique, 

Universite Laval. We thank the Canada Research Chair 

itt i r> u if • i- n ■ program and NSERC for support. 

We thank Roberto Maiolmo, Benjamin Oppenheimer, 

and Matthew Pieri for useful correspondence. The sub- 

APPENDIX 
THE METALLICITY OF THE GALAXIES 

Equation ([TJ provides the oxygen abundance versus mass and redshift. To calculate a total metallicity Zq, we need 
to make assumptions about the abundances of the other elements. The metallicity Zq is given by 



n z^z 

Z>2 



Zg = — ^ , (Al) 



Z>2 

where n z and \i z are the number density and atomic mass for element Z , respectively, and the sum is over all metals, 
that is all elements excluding hydrogen and helium. We rewrite equation (|A1|) as 

— Vz 

Z G = Z -^-^ = , (A2) 

n H/% + n Hc/%c + n O / , — n Z 
Z>2 ° 

where n Q is the number density of oxygen, and \i z is the equivalent oxygen atomic mass, defined by 

Vz.o ^/Zt^^z- (A3) 

Z>2 ° 

We assume that, when oxygen is produced, other metals are produced at the same rate, so the ratios n z /n Q remain 
constant. Of course, the actual ratios will depend on the relative importance of Type la and Type II supernovae, 
among other things, but this is a convenient approximation. We calculated fi z using the solar-system abundances 

n z given bv lCoxl (|2000f ). We only consider metals having an abundance log n z > 5 (where for oxygen, log n Q — 8.93), 
which includes the elements: C, N, O, Ne, Na, Mg, Al, Si, P, S, CI, Ar, K, Ca, Cr, Mn, Fe, and Ni. The resulting 
value is [i zo — 31.62380 
We then rewrite equation (|A2j) as 

n n zo n ot x z.o , A .v 

Z G = 7 r = : , (A4) 

„ (a 1 + n a n ^wM + n o^z,o 

n H I Mh + A*Hc I T n Of-Z,0 



n 



H 



where ji H is the equivalent hydrogen atomic mass, defined by 



AV.h = + — A*Hc ■ (A5) 
n„ 



We use the primordial abundances of H and He to calculate ^He/ n H- This is an approximation, since stellar evolution 
will increase the helium content relative to hydrogen. We get /« pr H = 1.3990. We can then rewrite equation (|A4[) as 

Zg = - ( . ^° • (A6) 

(n H /n )Mpr,H + Mz,o 

Using the ratio n H /n Q given by equation (JlJ, we can then calculate Zg for a given galaxy stellar mass M*. 



The main contributors are oxygen, iron, and carbon. 



18 



REFERENCES 



Aguirre, A., Hernquist, L., Schaye, J., Weinberg, D. H., Katz, N., 

& Gardner, J. 2001, ApJ, 560, 599 
Anders, E., & Grevesse, N. 1989, GeCoA, 53, 197 
Baldwin, J. A., Hamann, F., Korista, K. T., Ferland, G. J., 

Dietrich, M., & Warner, C. 2003, ApJ, 583, 649 
Barai, P., Gopal-Krishna, Ostcrman, M. A., & Wiita, P. J. 2004, 

BASI, 32, 385 
Barai, P., & Wiita, P. J. 2007, ApJ, 658, 217 
Barai, P. 2008, ApJ, 682, L17 

Bergeron, J., Aracil, B., Petitjean, P., & Pichon, C. 2002, A&A, 
396, Lll 

Bertone, S., Stoehr, F., & White, S. D. M. 2005, MNRAS, 359, 
1201 

Buries, S., & Tytler, D. 1996, ApJ, 460, 584 
Carr, B. J., Bond, J. R., & Arnett, W. D. 1984, ApJ, 277, 445 
Carswell, B., Schaye, J., & Kim, T.-S. 2002, ApJ, 578, 43 
Cen, R., & Bryan, G. L. 2001, ApJ, 546, L81 
Cen, R., & Ostriker, J. P. 1999, ApJ, 519, L109 
Chartas, G. et al. 2009, New Astronomy Reviews, 53, 128 
Cowie, L. L., Songaila, A., Kim, T.-S., & Hu, E. M. 1995, AJ, 
109, 1522 

Cox, A. N. 2000, Allen's Astrophysical Quantities (4th Ed.), AIP 

press, Springer- Verlag 
Crenshaw, D. M., Kraemcr, S. B., & George, I. M. 2003, 

ARA&A, 41, 117 
Danforth, C. W., & Shull, J. M. 2005, ApJ, 624, 555 
Dietrich, M., Hamann, F., Shields, J. C, Constantin, A., Heidt, 

J., Jagcr, K., Vestergaard, M., & Wagner, S. J. 2003, ApJ, 589, 

722 

Dietrich, M., Mathur, S., Grupe, D., & Komossa, S. 2009, ApJ, 
696, 1998 

Di Matteo, T., Croft, R. A. C, Springel, V., & Hernquist, L. 

2004, ApJ, 610, 80 
D'Odorico, V., Cristiani, S., Romano, D., Granato, G. L., & 

Danese, L. 2004, MNRAS, 351, 976 
Fabjan, D., Borgani, S., Tornatore, L., Saro, A., Murante, G. & 

Dolag, K. 2010, MNRAS, 401, 1670 
Ferrara, A., Pettini, M., & Shchekinov, Y. 2000, MNRAS, 319, 

539 

Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Furlanetto, S. R., & Loeb, A. 2001, ApJ, 556, 619 
Ganguly, R., & Brotherton, M. S. 2008, ApJ, 672, 102 
Germain, J., Barai, P., & Martel, H. 2009, ApJ, 704, 1002 
(Paper I) 

Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 
Gnedin, N. Y. 1998, MNRAS, 294, 407 
Gopal-Krishna & Wiita, P. J. 2001, ApJ, 560, L115 
Gopal-Krishna, Wiita, P. J., & Barai, P. 2004, JKAS, 37, 517 
Groves, B. A., Heckman, T. M., & Kauffmann, G. 2006, MNRAS, 
371, 1559 

Haiman, Z., & Loeb, A. 1997, ApJ, 483, 21 
Hamann, F., & Ferland, G. 1992, ApJ, 391, L53 
Hamann, F., Warner, C, Dietrich, M., & Ferland, G. 2007, 
ASPC, 373, 653 

Hcllsten, U., Dave, R., Hernquist, L., Weinberg, D. H., & Katz, 

N. 1997, ApJ, 487, 482 
Hinshaw, G., et al. 2008, ApJS, 180, 225 

Hockncy, R. W., & Eastwood, J. W. 1988, Computer Simulation 

Using Particles (New York: Adam Hilger) 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., 

Robertson, B., & Springel, V. 2006, ApJS, 163, 1 
Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 

731 

Juarez, Y., Maiolino, R., Mujica, R., Pedani, M., Marinoni, S., 

Nagao, T., Marconi, A., & Oliva, E. 2009, A&A, 494, L25 
Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 311, 576 
Khalatyan, A., Cattaneo, A., Schramm, M., Gottlober, S., 

Steinmetz, M., & Wisotzki, L. 2008, MNRAS, 387, 13 
King, A. R., & Pounds, K. A. 2003, MNRAS, 345, 657 
Kirkpatrick, C. C, Gitti, M., Cavagnolo, K. W., McNamara, B. 

R., David, L. P., Nulsen, P. E. J., & Wise, M. W. 2009, ApJ, 

707, L69 



Lapi, A., Shankar, F., Mao, J., Granato, G. L., Silva, L., De 

Zotti, G. & Danese, L. 2006, ApJ, 650, 42 
Levine, R., & Gnedin, N. Y. 2005, ApJ, 632, 727 
Lynden-Bell, D. 1969, Nature, 223, 690 
Madau, P., Ferrara, A., & Rees, M. J. 2001, ApJ, 555, 92 
Maiolino, R. et al. 2008, A&A, 488, 463 

Malbon, R. K., Baugh, C. M., Frenk, C. S., & Lacey, C. G. 2007, 

MNRAS, 382, 1394 
Mathur, S., & Fields, D. 2009, in Future Directions in Ultraviolet 

Spectroscopy, AIP Conference Proceedings, Vol. 1135, p. 209 
Matsuoka, K., Nagao, T., Maiolino, R., Marconi, A., & Taniguchi, 

Y. 2009, A&A, 503, 721 
Mcnci, N., Fiore, F., Puccetti, S., & Cavaliere, A. 2008, ApJ, 686, 

219 

Moll, R. et al. 2007, A&A, 463, 513 

Monaco, P., & Fontanot, F. 2005, MNRAS, 359, 283 

Nagao, T., Marconi, A., & Maiolino, R. 2006a, A&A, 447, 157 

Nagao, T., Maiolino, R., & Marconi, A. 2006b, A&A, 447, 863 

Nesvadba, N. P. H., Lehnert, M. D., Eisenhaucr, F., Gilbert, A., 

Tecza, M., & Abuter, R. 2006, ApJ, 650, 693 
Nesvadba, N. P. H., & Lehnert, M. D. 2007, Heating versus 

Cooling in Galaxies and Clusters of Galaxies, ESO Astrophysics 

Symposia, (Springer- Verlag Berlin Heidelberg), 416 
Nesvadba, N. P. H., Lehnert, M. D., De Breuck, C, Gilbert, A. 

M., & van Breugel, W. 2008, A&A, 491, 407 
Oppenheimer, B. D., & Dave, R. 2006, MNRAS, 373, 1265 
Oppenheimer, B. D., Dave, R., & Finlator, K. 2009, MNRAS, 

396, 729 

Ostriker, J. P., & Gnedin, N. Y. 1996, ApJ, 472, L63 

Pieri, M., Martel, H., & Grcnon, C. 2007, ApJ, 658, 36 

Prescott, M. K. M., Dey, A., & Jannuzi, B. T. 2009, ApJ, 702, 554 

Rauch, M., Haehnelt, M. G., & Steinmetz, M. 1997, ApJ, 481, 601 

Rees, M. J. 1984, ARA&A, 22, 471 

Richstonc, D. et al. 1998, Nature, 395, A14 

Ryan-Weber, E. V., Pettini, M., Madau, P., & Zych, B. J. 2009, 

MNRAS, 395, 1476 
Salpeter, E. E. 1964, ApJ, 140, 796 

Sazonov, S. Y., Ostriker, J. P., Ciotti, L., & Sunyaev, R. A. 2005, 

MNRAS, 358, 168 
Scannapieco, E., & Broadhurst, T. 2001, ApJ, 549, 28 
Scannapieco, E., Ferrara, A., & Madau, P. 2002, ApJ, 574, 590 
Scannapieco, E., & Oh, S. P. 2004, ApJ, 608, 62 
Scannapieco, E., Pichon, C, Aracil, B., Petitjean, P., Thacker, R. 

J., Pogosyan, D., Bergeron, J., & Couchman, H. M. P. 2006, 

MNRAS, 365, 615 
Schaye, J., Rauch, M., Sargent, W. L. W., & Kim, T.-S. 2000, 

ApJ, 541, LI 

Schaye, J., Aguirre, A., Kim, T.-S., Theuns, T., Rauch, M., & 

Sargent, W. L. W. 2003, ApJ, 596, 768 
Shankar, F. 2009, New Astronomy Reviews, 53, 57 
Shen, S., Wadsley, J., & Stinson, G. 2009, MNRAS, submitted 

(arXiv: 0910.5956) 
Shull, J. M., Tumlinson, J., & Giroux, M. L. 2003, ApJ, 594, L107 
Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 2002, ApJ, 578, 

737 

Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 2004, ApJ, 606, 
92 

Songaila, A., & Cowie, L. L. 1996, AJ, 112, 335 
Songaila, A. 1997, ApJ, 490, LI 
Songaila, A. 2001, ApJ, 561, L153 

Tegmark, M., Silk, J., & Evrard, A. 1993, ApJ, 417, 54 

Thacker, R. J., Scannapieco, E., & Davis, M. 2002, ApJ, 581, 836 

Tornatore, L., Borgani, S., Viel, M., & Springel, V. 2010, 

MNRAS, 402, 1911 
Tripp, T. M. et al. 2002, ApJ, 575, 697 
Voit, G. M. 1996, ApJ, 465, 548 

Wang, T., Zhou, H., Yuan, W., Lu, H., Dong, X., & Shan, H. 

2009, ApJ, 702, 851 
Wise, J. H., & Abel, T. 2008, ApJ, 685, 40 



